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ABSTRACT 

An  extension  of  Wetzel' s  theory  is  proposed  to  explain  electron 
precursor  phenomena  in  pressure  driven  shock  tubes.  The  theory  is 
based  on  the  assumption,  which  has  been  experimentally  verified,  that 
radiation  from  the  hot  gas  behind  the  shock  front  is  responsible  for  the 
precursor  effect.  A  point  radiator,  located  at  the  shock  front  and 
moving  with  the  shock  velocity,  is  assumed.  It  is  shown  that  this 
radiator  can  be  idealized  by  a  black  body.  The  precursor  electron 
density  is  derived  from  a  one -step  photoionization  of  the  impurities 
present  in  the  driven  gas. 

As  an  application  of  the  theory,  precursor  profiles  are  calculated 
for  an  impurity.  The  partial  pressure  of  i*  •  varied  from  2  to  75^iHg. 
The  temperature  behind  the  shock  is  assumed  to  be  11,  000°K. 
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SECTION  I 


INTRODUCTION 

Upon  reentry  of  a  satellite  or  meteor  into  the  earth' s  atmosphere, 
many  investigators  observed  ionization  far  in  front  of  the  shock  wave 
of  the  reentering  body.  S.  C.  Lin*  and  his  co-workers,  through  their 
radar  measurements,  established  that  high  electron  densities  are  found 

2 

several  meters  ahead  of  the  shock  of  a  manned  satellite.  Cook  and  Hawkins 
also  found  detectable  electron  densities  in  front  of  reentering  meteors. 

Prior  to  these  observations,  a  number  of  investigators  working  with  shock 
tubes  observed  signals  on  probes  far  in  front  of  the  shock.  All  these  signals 
became  known  as  precursor  phenomena. 

The  precursor  ionization  is  important  for  a  number  of  reasons. 

First,  the  problem  of  discrimination  of  reentering  bodies  may  be  highly 
affected  because  this  phenomena  may  change  measurably  the  radar  cross- 
section  of  reentering  bodies.  Secondly,  since  the  precursor  ionization 
must  be  associated  with  energy  transfer  from  behind  the  shock  to  in  front 
of  the  shock,  the  conventionally  assumed  shock  structure  may  have  to  be 
modified  to  take  into  account  this  energy  transfer.  Thirdly,  precursor 
ionization  causes  the  shock  .  >nt  to  pass  through  a  partially  ionized  gas, 
and  this  may  affect  the  chemical  reactions  taking  place  behind  the  shock. 

Several  hypotheses  have  been  advanced  to  explain  precursor 
phenomena.  According  to  one  hypothesis,  the  shock  heated  gas  radiates 
U.  V.  and  X-rays  which  photoionizes  the  gas  in  front  of  the  shock^.  This 
ionization  may  be  a  one -step  process  from  the  ground  state  of  the  gas, 
or  it  may  be  a  two-step  process  where  the  atom  is  first  excited  and  then 
ionized.  Also  the  source,  the  shock  heated  gas,  may  behave  as  a  continuum 
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black  body  radiator  or  line  radiation  may  be  significant.  Recently, 

4  5 

Clarke  and  Murty  have  suggested  that  line  radiation  plays  a  significant 
role  in  precursor  phenomena. 

Another  hypothesis  begins  with  the  recognition  that  a  large  electron 
and  ion  density  gradient  exists  at  the  shock  front.  This  would  result  in 
a  strong  tendency  for  the  electrons  and  ions  to  diffuse  into  the  precursor 
region.  Weyman^  claimed  experimental  verification  for  electrons  escaping 
from  behind  the  shock  with  velocities  several  times  larger  than  the  shock 
velocity.  A  critical  theoretical  study  of  this  diffusion  hypothesis  by 
Petschek^  and  Wetzel^  concluded  that  this  phenomena  can  be  completely 
neglected  in  the  precursor  region  except  in  a  region  about  one  Debye 
length  ahead  of  shock  front.  The  theoretical  models  mentioned  above 
do  not  involve  the  gas  dynamics  of  the  shock  in  any  essential  way*  i.  e. , 
the  basic  equations  of  fluid  mechanics  are  uncoupled  from  the  equations  of 
radiation  transport. 

Besides  these  two  hypotheses,  there  are  some  other  postulations 

8  9 

suggested  to  explain  this  ionization  phenomena.  Hollyer  and  Gloersen 
attribute  the  precursor  to  the  effect  of  photoemission  from  the  walls 
of  shock  tube.  Photoemission  would  leave  an  insulated  shock  tube  wall 
positively  charged  and  only  electrons  would  be  collected  by  electrostatic 

probes.  However,  the  presence  of  la:-ge  numbers  of  positive  ions  have 

,  .  .  .  .10 

been  found  in  the  precursor  region  . 

Lederman  and  Wilson^  applied  a  microwave  resonant  cavity 
technique  to  measure  precursor  profiles.  This  technique  avoids  the 
difficulties  and  uncertainties  of  interpreting  probe  signals.  The  cavity 
technique  yields  definite  values  of  electron  density  while  the  probe  measure¬ 
ments  give  only  the  shape  of  the  precursor  profile  and  not  always  the  absolute 
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electron  number  density.  Lederman  and  Wilson' s  experiments  concluded 
that:  (1)  the  precursor  ion’zation  is  due  to  radiation  from  the  shocked  gas 
behind  the  incident  snock;  (2)  the  electron  precursor  velocity  is  equal  to 
the  shock  velocity. 

1 1  1 

Hammerling  and  S.  C»  Lin  analyzed  the  precursor  phenomena 
for  a  satellite  reentering  the  atmosphere  by  assuming  that  the  ultra¬ 
violet  radiation  comes  from  the  b,l£  +-«  X1  transition  of  N£*  This 

u  g 

radiation  then  precedes  the  shock  into  the  unshocked  gas  and  ionizes  the 
molecular  oxygen.  They  found  that  the  precursor  ionization  halo  appears 
satisfactory  for  explaining  the  leading  edge  echo  before  the  MA-6  Mercury 
capsule.  It  is  worthy  to  note  that  their  model  is  based  on  a  one  step  ionization 
process. 

So  far,  the  important  role  that  impurities  may  play  in  precursor 

ionization  has  not  been  emphasized  in  the  literature.  This  is  true  even 
7  12 

though  Petschek  and  Sturtevant  recognized  that  impurities  play  a 
dominating  role  in  initiating  ionization  in  high  temperature  gases  behind  shocks. 

In  this  paper  a  one-step  photoionization  theory  will  be  presented. 

This  theory  is  based  on  a  Wetzel -type^  model  with  several  important 
modifications.  The  important  role  that  impurities  play  in  precursor 
ionization  will  be  emphasized. 

The  authors  would  like  to  thank  Prof.  Samuel  Lederman  for  many 
enlightening  discussions  during  the  course  of  this  work. 
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SECTION  II 


PHOTOTONIZA1  ION  THEORY  OF  PRECURSOR  IONIZATION 

The  mechanism  proposed  by  Wetze’  '  views  the  shock  wave  as  a 
radiator  moving  at  the  shock  speed;  this  radiator  ionizes  the  gas  ahead 
of  the  shock,  producing  the  precursor.  The  experiments  of  Ref  10 
suggest  this  type  of  mechanism. 

II.  A.  Review  of  Wetzel*  s  Theory 

A  one -dimensional,  shock-based  coordinate  system  is  assumed 
in  Fig.  1.  In  this  coordinate  system  the  gas  is  streaming  in  the  minus 
x-direction  at  the  shock  velocity,  U.  Steady  state  conditions  are  assumed. 
The  governing  equations,  momentum  trans.. ort  equations  for  electrons  and 
for  ions  and  Poisson' s  equation  to  provide  electrostatic  coupling  between  the 
two  charged  species,  follow: 


D 

e 


D. 

1 


d2n 


dx 


e  cx 


d2n 


dx 


1  -  u  d 

1  ^  “cDT 


an 

(En  )  +  U  — —  =  -  r(x)  +  s(x) 
®  d.x 


dn. 

(En. )  +  U  - —  =  -  r(x)  +  s(x) 

1  dx 


(1) 

(2) 


co  “3T  =  e(ni  '  (3) 

where  n  and  n.  are  the  electron  and  ion  densities,  and  D  ,  U  an  1  D.,  u.  are 
e  i  e  e  i  i 

the  respective  diffusion  and  mobility  coefficients  of  electrons  and  ions  relative 
to  the  flow  velocity,  U.  E  is  the  electrostatic  field,  and  r(x)  and  s(x)  are 
source  and  sink  functions,  respectively.  Eq.  (1),  for  example,  is 
simply  a  balance  equation  for  electrons.  The  first  term  on  the  left  is  a 
net  loss  due  to  diffusion,  the  second  term  a  loss  due  to  electric  forces,  and 
the  third  is  the  familiar  conduction  term;  this  is  balanced  on  the  right  by 
sink  and  source  terms. 
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The  sink  term,  s(x),  is  due  to  losses  such  as  recombination.  Such 
losses  are  unimportant  in  the  present  application,  and  s(x)  will  therefore 
be  dropped  from  the  equations.  The  source  term,  r(x),  introduces 
additional  electrons  and  ions  due  to  photoionization,  and  it  is  this  term 
which  has  to  be  further  investigated. 

We  define  a  partial  source  function,  r(x;v)  dv.  as  the  number  of  ioni¬ 
zations  per  unit  volume  per  second  produced  at  x  by  photons  in  the  frequency 
range  from  v  kj  v  +  d\>  With  each  r(x;v)  dv  we  will  associate  a  partial 
ionization  density,  n(x;v)dv,  with  each  range  of  photon  frequencies.  The 
total  source  and  density  functions  then  become 

oo 

(• 

r(x)  =  J  r(x;v)dv 
vi 

(4) 

00 

* 

n(x)  =  \  n(x;v)dv 
“vi 

where  m.  is  the  lowest  ionization  frequency  of  the  gas  in  the  precursor  region. 

Now  let  N(v)  be  the  number  of  photon*  emitted  per  unit  area  per 
second  per  unit  frequency  range  by  the  radiator  located  at  x  =  0.  The  prob¬ 
ability  that  a  photon  of  frequency  v  will  penetrate  a  distance  x  and  be  ab¬ 
sorbed  between  x  and  x  +  dx  is  then 

ng  Q(v)  exp  J-  ng  Q(v)  x  ]dx  (5) 

where  n  is  the  number  density  of  the  gas  particles  in  the  precursor  region 
and  Q(v)  is  the  absorption  cross-section  of  photons  of  frequency  v  in  the  gas. 

If  Q^(v)  is  the  ionization  cross-section  for  photons  of  frequency  v»  then 
Q^(v)/Q(v)  is  the  probability  that  an  absorbed  photon  will  produce  an  ioniza¬ 
tion,  and  r{x;v)  is  given  by 
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(6) 


r(x  ;v)  =  N(v)  ng  Q.(v)  exp 


[•  "6q<v)  x] 


From  the  experimental  results  in  Ref.  10,  we  can  neglect  all  terms 
involving  since  the  precursor  profile  is  caused  by  photoionizat'on 

only.  Thus,  the  first  two  termj  in  Eqs.  (1)  and  (£)  will  be  dropped, 
leaving  simply 

.  U  dn-fo-V.)  =  N(v)  n  Q  (v)  exp  [-  n  Q(v)  xl  =  r(x;v)  (7) 

dx  *  L  o  J 


where  a  distinction  betwen  electron  and  ion  distribution  is  no  longer  neces¬ 
sary.  Now,  we  will  assume  that  the  shock  is  moving  at  nearly  constant 
velocity.  This  implies  that  the  attenuation  of  shock  strength  is  negligible. 
This  assumption  is  also  verified  experimentally. 

We  can  integrate  Eq.  (7)  with  respect  t  x  from  oo  to  x  using  the 
boundary  condition  that  n(oo;v)  =  0;  and  also  integrating  with  respect  to  v, 
we  find  the  total  ionization  density: 

.1  -00  Q4(v) 

n(x)  =  U  j  N(v)  -  exp 

J  Vi  Q(v) 

To  further  simplify  Eq.  (8),  let  us  define  f(v)  =  Q(v)/Q(\i)  as  the 
absorption  cross-section  normalized  to  its  value  ^t  the  ionization  edge,  and 

X 

7}  ~~~  as  a  dimensionless  distance  upstream  of  the  radiator  (where 
1 

-  1/n^  Q(v  )  is  the  photon  mean  free  path  at  the  ionization  edge).  In  terms 
i  g  i 

of  f(v)  and  r] ,  Eq.  (8)  can  be  written 

oo 

n(rj)  =  U  1  J  N(v)  f(v)  1  exp  \-  f(v)n  ]  dv  (9) 


[-r.gQ(v)x  ], 


(8) 
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II.  B.  Precursor  Distribution  in  Argon 


It  is  first  assumed  that  the  shock  frc  nt  radiates  as  a  black  body  at 
a  temperature  equal  to  the  temperature  behind  the  sho-ik,  ar.d  that  it  is  the 
high  energy  photons  of  this  radiation  that  are  responsible  for  the  precursor. 
Using  the  high  energy  limit  of  the  black  body  curve,  we  get 


N(v)  = 


Zi 


2 

v  exp 


hv 

"KT 


,  hv  »  KT 


(10) 


For  many  absorption  processes,  it  can  be  assumed  that  each  absorbed 
photon  produces  an  ion  pair.  For  such  processes  we  have  Q.(v)/Q(v)  =  1» 

o 

This  is  approximately  correct  for  inert  gases,  such  as  Argon.  From  790A 

o 

(which  corresponds  to  the  first  ionization  potential  of  Ar ?)  to  about  420A, 

the  cross-section,  Q(v),  for  Argon  is  fairly  constant  and  equal  to  about 

- 18  2  0 
35  x  10  cm  .  For  wave  lengths  smaller  than  420A,  the  cross-section 


falls  off  approximately  as  v 
values  of  f(v)  are  used: 

f(v)  =  1 

f(v)  = 


-3 


Therefore,  in  the  calculation,  the  following 


for  420A  <  X 


790A 


V420A 


(ID 


for 


X  <  420A 


Using  Eqs.  (10)  and  (11)  inEq.  (9),  one  cnn  find  n(rj)  as  a  function 
of  q  .  The  calculated  distribution  falls  off  much  too  rapidly  to  explain  the 
observed  precursor  *  For  example,  at  50  cm.  ahead  of  shock,  the 

calculated  number  density  is  5  orders  of  magnitude  below  the  observed 

,  10,  12 

value 

The  electron  number  density  at  the  shock  front  is  obtained  by  letting 
r|  -*  0  in  Eq.  (9): 
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SECTION  III 


MODIFICATIONS  IN  THE  PHOTOIONIZATION  THEORY 

Wetzel,  as  described  in  Section  IIA,  assumes  a  plane  source,  infinite 
in  extent,  in  deriving  his  relations  for  the  precursor  profile  ahead  of  the 
shock.  We  do  not  feel  this  is  a  realistic  assumption  to  make  for  any  practical 
confirmation  of  his  theory,  the  source  is  seldom  of  this  type.  For  example, 
for  a  reentering  vehicle,  the  source  behaves  as  a  point  radiator  located 
at  the  nose  of  the  vehicle,  and  some  account  must  be  taken  of  the  fact  that 
the  radiation  falls  off  as  the  inverse  aqua  re  of  the  distance  from  the  vehicle 
even  if  there  is  no  absorption  of  the  radiation.  Also,  in  the  shock  tube,  the 
diameter  of  the  shock  tube  is  small  compared  to  the  distance  ahead  of  the 
shock  where  measurements  are  of  interest,  and  it  is  more  realistic  to 
assume  a  point  source  at  the  shock  front  than  a  plane  radiator. 

In  thi;;  paper  a  point  radiator,  located  at  the  shock  front,  is  assumed. 

At  the  frequencies  of  interest,  it  can  rafely  be  assumed  that  all  radiation 
impinging  on  walls  of  the  shock  tube  is  completely  absorbed.  Since  the 
radiation  will  now"fall  off"  faster  than  for  a  one -dimensional  radiator,  a 
one-step  photoionization  of  Argon  will  still  not  supply  the  answer  to  what 
is  responsible  for  the  precursor.  this  paper,  the  important  role  that 
impurities  in  the  driven  gjs  of  the  shock  tube  have  on  the  nature  of  the 
precursor  will  be  emphasized. 

III. A.  The  Electron  Precursor  Distribution 
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that  our  m-a.ur.ment.  are  made  far  enough  from  the  .ource  that  the 
radiation  flux  is  everywhere  in  the  axial  direction. 

Thi.  is  inserted  into  Eq.  (7),  and  the  integration  is  again  from  to  to  x: 


n(x:s)  =  -  — n  N(v)  Qi<v)  J 


exp 


4nU 


OO  X 


-  n 

£ 


Q(v)x 


dx 


Integrating  Eq.  (14)  by  parts,  we  get 
n(x;v)  =  -  - — —  n^N(v)iQi(v) 


4tt  U 


h  exp  [•  “( 

. . L,  4“p  [-“g 


Q(v)  * 


dx 


The  integral  on  the  right  side  of  Eq.  (15)  can  be  written  as  follows 


-n  Q(v)x 


g 

-oo 


^  i  «p  [-  “g  Q(v)  *  H  ■  J  _ 

=  Ei  [-  ng  Q(v)]x  =  Ei(  -  ?) 


(14) 


(15) 


( 16> 


where  ?  =  n  Q(v)  *• 
g 


14,  15. 


The  function  Ei(-  5)  is  the  well-known  exponential  integral  ;its 

series  expansions  can  be  written  as: 


co 


n.n 


Ei(-§)  =  y  +  In  §  +  £  (  n(n!  j  *  (?  >Q) 

n*l 


(17) 


wnere 


v  =  iim 

n-»  0 


>+7  ♦  T  +  *  *  *  *  +  5  -l08n 


=  0.577 


(y  is  called  Euler  or  Mascheroni' s  constant) 
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For  5  >  >  1,  El(  -  {)  has  the  following  asymptotic  representation: 


Ei(  -  §)  = 


c-5 

T?) 


1*  2* 
1  -  —  +  "i 


?  s 

In  terms  of  Ei(  -  5),  E3q.  (15)  can  be  written 

A 


3'. 
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^  *  »  •  » 


n(x;v)  =  - 


4tt  U 


■n  N(v)^(v) 


] 

1 

-  —  exp 

-  ng  Q(v)  x 
.  ^ 

ng  Q(v)  Ei(-  5)[ 


(18) 


(19) 


The  total  number  density,  n(x),  is  obtained  by  integrating  Eq.  (19) 
with  respect  to  v« 


A  n 

n(x)  =  - S_ 

4tt  U 


J 


•°°  N(v) 


exp 


00 


ng  Q(v)  x 


dv 


+  n  J  N(v)Qi(v)Q(v)  Ei(  -  5)  dv 


(20) 


Equation  (20)  now  becomes  the  revised  form  of  Eq.  (9)  in  the  new 
theory.  It  should  be  mentioned  again  that  Eq.  (20)  will  cause  an  even 
faster  decrease  with  distance  in  the  precursor  electron  concentration 
than  given  by  the  theory  of  Section  II A,  since  the  source  radiation  falls 
off  faster.  For  this  reason,  the  conclusions  reached  in  Section  IIB  con¬ 
cerning  an  Argon  photoionization  interpretation  of  the  precursor  phenomenon 
will  not  change. 

III.B.  The  Black  Body  Assumpticn 

The  radiant  energy  required  for  photoionization  of  the  gases  ahead 
o*  hock  waves  comes  from  theshock  heated  high-temperature  gas  behind 
the  shock  front.  Thus,  it  is  necessary  to  know  the  time  history  of  the  light 
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intensity,  the  spectrum,  and  the  absorption  coefficients  of  radiating  sources 
before  an  accurate  prediction  of  the  precursor  profile  is  possible.  It  is 
difficult  to  analyze  this  non-steady,  non- equilibrium  radiating  source  behind 
the  shock  analytically  and, unfortunately,  spectrum  data  for  shock  heated  gases 

are  also  virtually  non-existent. 

7  16 

Petsfhek  and  Byron  and  Blackman  and  Niblett  studied  relaxation 
processes  behind  Argon  shocks  and  found,  experimentally,  that  the  relaxation 

time  was  strongly  dependent  on  the  impurity  content  of  the  Argon  driven  gas. 

12 

Recently,  Sturtevant  applied  a  mass  spectrometer  to  study  the  ionization 
processes  and  chemical  species  behind  an  Argon  shock.  He  found  a  large 
amount  of  0+  and  H*  (about  five  times  more  abundant  than  Argon)  during 
the  initial  stages  of  ionization.  The  6mall  amount  of  impurities,  he  concluded, 
played  a  dominant  role  in  the  ionization  process,  or  at  least  in  the  initial 
ionization. 

In  Section  IIAthe  precursor  ionization  profile  is  calculated  by  assuming 
that  the  shock  is  a  black  body  radiator,  i.  e. ,  the  radiating  source  is  optically 
thick.  This  is  a  reasonable  assumption  to  make  for  frequencies  above  the  first 
ionization  edge  of  Argon.  However,  for  frequencies  below  this  ionization  edge, 
the  radiation  is  not  necessarily  expected  to  behave  like  black  body  radiation. 

Due  to  impurities  in  the  driven  gas,  it  may  still  be  safe  to  make  a  black  body 
as  mption.  An  attempt  will  be  made  in  the  present  section  to  answer  this 
question. 

The  equation  for  radiative  transfer  in  steady  case  at  the  absence  of 
scattering  is 


-hi) 
V  V 


(21) 
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where 

pjv  =  emission  coefficient 
k^  =  absorption  coefficient 

I  =  intensity  of  radiation 

v 

Under  the  assumption  of  thermodynamic  equilibrium  in  the  shocked 
heated  gas,  Bq.  (21)  reduces  to 


3  I 

- =  k  (B  -  I  ) 

az  v  v  v 


(22) 


where 

B  (v)  = 

V  >/kT-l 

After  integrating  Eq.  (22)  and  using  the  boundary  condition  that 

I  =  0  at  Z  =  0  (at  the  contact  surface),  we  get 
v 

r  -k  Z, 

'v  =Bv  [*•■  V  ] 


(23) 


(24) 


The  total  radiation  intensity  is 


I  =  (*  B  (v)  (1  -  e  v  )  dv  (25) 

ihe  unknown  in  Eq.  (25)  is  the  absorption  coefficient  "k  "  which 
is  a  function  of  temperature  as  well  as  frequency.  If  k^  Z  »  1,  it  is 
optically  thick  and  1^  will  be  very  similar  toa  black  body  emitter. 

There  are  a  variety  of  mechanisms  for  emission  and  absorption 
of  radiation  which  are  associated  with  transitions  between  energy  levels 
of  the  gas  atoms  or  molecules.  In  these  transitions,  the  change  in  in¬ 
ternal  energy  of  the  atom  is  equal  to  the  radiant  energy  absorbed.  The 
types  of  transitions  that  are  important  in  a  given  situation  depends  on  the 
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nature  of  the  gas,  on  the  distribution  of  molecules  among  the  various  energy 
levels,  and  on  the  radiant  frequency  range  of  interest.  They  cai  be  roughly 
classified  as  3  kinds,  i.  e.;  bound  -bound,  bound-free,  and  free  free  transitions. 
Bound- bound  transition  gives  rise  to  a  discrete  line  spectrum  (  the  contri¬ 
bution  of  this  type  of  transition  will  be  neglected  in  this  report),  while  free- 
free  and  bound-free  transitions  give  rise  to  a  continuous  spectrum. 

The  composite  gas  molecules  in  the  high  temperature  region  consists 
mainly  of  Argon  with  some  impurities  or  trace  gases  such  as  dissocial. :d  air 
molecules  (oxygen  or  nitrogen),  etc.  However,  the  exact  amount  of  impuri¬ 
ties  mixing  with  the  Argon  is  not  known. 

Under  qua  si -equilibrium  conditions,  the  absorption  coefficient  k  can 
be  estimated  by  the  existing  theories.  The  first  step  is  to  evaluate  the 
population  of  various  atomic  and  molecular  energy  levels.  This  includes 
the  calculation  of  dissociation  and  the  degree  of  ionization.  The  second, 
a  problem  in  quantum  theory  of  radiation,  is  the  calculation  of  the  radiative 
transition  probability  for  transition  between  the  different  energy  levels.  This 
problem  is  very  difficult,  and  exper' mental  data  will  be  used. 

The  primary  impurity  in  the  driven  gas  is  assumed  to  be  (15yiHg). 

For  an  initial  pressure  of  Pj  =  2  mm  Hg  Argon  and  a  Mach  number  of  13.  2, 
the  temperature  behind  the  incident  shock  \  ,-proximately  11,  000°K.  The 
density  jump  acrors  the  shock  is  assumed  to  be  the  same  as  for  the  Argon, 

and  the  behind  the  incident  shock  is  assumed  to  be  fully  dissociated. 

17 

The  table  gives  the  calculated  values  for  k  for  atomic  oxygen  for  wave 

00  v 

lengths  between  930A  and  600A. 
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Table:  Absorption  Coefficients  of  Atomic  Oxygen 
Pj  (argon)  =  2  mm  Hg. 

M  -  13. 2 


Initial  O.,  Pressure  =  15|iHg. 


\(A) 

930 

900 

mm 

830 

800 

770 

730 

700 

600 

k  (cm  *) 

V 

4.  1  x  10"2 

4.  lx  10'2 

4.  2x10 ’2 

4. 6x10  ^ 

4.  9xl0'2 

8.  6xl0_1 

1.  lxlO'1 

3BSS 

A  typical  value  of  k^  for  atomic  oxygen  is  4  x  10  2  cm  For  a 
radiating  source  of  10  cm  in  length  (distance  between  shock  front  and  contact 

surface),  k  x  =  .  4. 

v 

Thus,  from  Eq.  (24),  we  get: 

=  VI  -e°-4)  =  0.35Bv 

Ihus,  the  photon  emitter  is  about  35%  efficient  and  can  be  idealized  as  a 
black  body  emitter. 
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III.  G.  Variation  in  the  Precursor  Profile  with  Impurity  Partial  Pressure 

The  theory  predicts  that  a  plot  of  precursor  electron  number  density 
versus  impurity  partial  pressure,  for  fixed  distance  ahead  of  shock,  should 
peak  at  a  certain  value  of  impurity  partial  pressure.  This  can  be  seen  from 
the  following  analysis: 

If  a  maximum  value  of  n(x)  occurs  at  a  certain  impurity  partial 
pressure,  p,  then 


and 


±^<o.r-h 

3p 


(26a) 


(26b) 


where  p^  =  the  partial  pressure  of  impurity  which  will  give  the  ionization 
profile  a  peak  value. 

We  differentiate  Eq*  (20)  with  respect  to  p,  and  recognize  that 
Q(v),  <^(v).  v  and  x  are  independent  of  p,  while  ng  is  not.  The  relation¬ 


ship  between  p  and  n  ist 

g 


p  =  n  KT. 
*  g 


When  we  set  the  derivative  equal  to  zero,  we  have 


(27) 


&  n(x)  _  3  n(x)  3  ng  _  1  3  n(x)  _ 

Tp - —  Tp  '  TTT  “Tn - ° 

g  g 

for  KT  ={=  0,  so  _  0 


(28) 

(29) 


get 


Performing  the  differentiation  in  Eq.  (20)  with  respect  to  n  ,  we 

g 


J_n (x) 

3  n 

g 


n 

g 


N(v)  Q(v)  Q.(v) 


+  2  Ei(-?) 


dv 


0 


(30) 
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The  second  derivation  of  n(x)  with  respect  to  "p"  follows: 


2  2  00 
9  n(x)  f  1  1  " 


=  f—  J  [N{V)Q(v)Q.(v)]  [  e'?  +  2Ei(-?)]  dv  (31) 

Sp  LktJ  v  1 


We  are  going  to  show  that - <  0  for  certain  values  of  n  or  p. 

&P  8 


The  exponential  integral,  Eif-^),  is  negative  in  value,  i.  e. ,  'i(  -  £)  <0; 

also  N(v)Q(v)Q.(v)  >  0,  ng>0and(kT)a  >  0. 

juppose  we  impose  0  <  5  <1  with  the  range  of  frequency  <  v  <  00; 

if  dn(x)/d  n  =  0  is  true,  then  comparison  of  Eqs,  (30)  and  (31)  leads  to  the 
8 

2  2 

result  that  9  n(x)/8p  <  0.  This  means  at  least  one  relative  maximum  value 

for  n(x)  occurs  somewhat  within  0  <5  <1.  On  the  other  hand,  if  f  >1  and 

2  2 

8n(x)/S  ng  =  0  are  true,  then  we  can  see  8  n(x)/&p  >  0  from  Eqs.  (30)  and 
(31)  ,  *0  that  at  least  a  relative  minimum  value  of  n(x)  occurs  somewhere 
at  ^  >  1.  In  what  follows,  it  will  be  seen  that  5  is  usually  restricted  to 
values  between  0  and  1  and,  therefore,  we  would  expect  to  find  maximum 
values  for  n(x). 
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PPk  jIONI zation  of  a  molecular  hydrogen  impurity 


IV.  A.  "Peak  Behavior"  for  a  Molecular  Hydrogen  Impurity 

To  evaluate  the  value  of  p  j-  n  that  will  give  the  precursor  profile 

S 

a  peak  value  (or  Sn(x)/3n  =  0),  we  have  to  solve  the  integral  Eq.  (30). 

S 

Tnis  equation  car.  be  solved  by  an  asymptotic  expansion. 

A  black  body  emitter,  i.  e.  ,  Eq.  (10)  will  be  used  and  the  impurity 
will  now  be  assumed  to  be  Molecular  Hydrogen,  whose  absorption  cross- 
nection  can  be  expressed  in  closed  form  as  follows: 

Q(v)  _  Vi  * 

TTT^)  "  “ 

where  hv.  =  15.422  ev. 


(32) 


Then  Eq.  (30)  can  be  written  as 

CD 


8n(x) 

Tn 

g 


2rr  f 

— T  Q'vi-  I 

c 


f(v)  —  exp 


i  i 


hV; 


KT 


v 

vi 


V 


'i  KT 


tz: 


) 


+  2n  Q(Vi)  ^  f(v)  Ei(-  e)  exp 


KT 


i  I  V  x 


dv  =  0 


(33) 


This  analysis  is  valid  for  any  impurity  whose  absorption  cross-section 
has  this  frequency  dependence.  When  Argon  is  the  driven  gas  and 
Molecular  Hydrogen  the  impurity,  corrc  ,tions  to  this  analysis  must  be 
made.  These  corrections  do  not  radically  change  the  calculated  profiles. 
This  will  be  fully  discussed  in  Ref.  18,  where  a  comparison  between 
theory  and  experiment  is  made. 
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where 


-1 


X.  =  C  n  Q(Vi)]  . 


Q4(v) 

f(v)  “STvT  ’ 


and 


§  =  K  Q(v)  X. 


The  saddle  point  of  integration  in  Eq.  (33)  is  at 


d  fv  ,  x  i  KT  I  _  „ 


cfP- 


f“ 

lvs 


\ .  V  tl  V-  J 
1  1 


(34) 


or 


v  ,  "““KT' 

"vT  =  V  XT  ffvT 


(35) 


Furthermore,  we  will  let  T  =  11,  00c  K  (  the  temperature  behind  the 
shock),  so  that 


hv4 

TT  =  16,3 

Also  since  x/\.  <  1,  v/v^  is  less  than  l  and  the  saddle  point  is  located 

outside  the  range  of  integration.  Thus  Eq.  (33)  can  be  expanded  asymptotically 

19 

by  using  the  formula  : 


00 


f  G(s)  e  nT{8)  ds  =  -  1 


G(s  )  e01*8^ 


n 


t  '  (8a) 


CD 

•  J_  I  d 

n  J  "3s 

sa 


G(s) 


t'(s) 


enT<8)  ds 


(36) 


where  Ois  a  numerical  parameter  which  determines  the  accuracy  of  the 
convergence  of  the  expansion.  In  our  case,  n  =  h  v./KT  =  16.  3.  After 


applying  formula  (36)  consecutively  two  times,  and  after  considerable 
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simplification,  we  get 


where 


1.  06 
*1 


-  0.  058  .  2  e^  Ei{  ~  n)  +  0{  — - - )  =  0 

lb-  33 


T)  =  ng  QtV^X. 


(37) 


Equation  (37)  has  been  solved  graphically  and  its  solution  is 
r  ~  0.  C.  .  The  absorption  cross-section  for  H_  at  the  ionization 

in&Xi  Lt 

edge  is 

Q(v.)  =  10‘17cm2  (38) 

Using  the  aoove  values,  the  partial  pressure  of  at  which  the  peak 
occurs  is  obtained  as  a  function  of  distance  ahead  of  the  shock  ana 
is  plotted  in  Fig.  2. 


IV.  B.  The  Precursor  Profile 

If  a  black  body  source  (Eq.  10)  ie  assumed,  integration  of  E~.(20) 
(repeated  use  of  Eq.  (36)is  made)gives  the  following  expression  for  the 
precursor  profile: 


n(x)  = 


3  - 16.  3  -  tj 

Av.  n  Q(v) 

_ 1 _ S _ i_  e 

2  U  c2x 


1 

1573 


[  1  +  e 


?! 


Ei(-q  )] 


+  ( 


1 

IE73 


) 


2 


(39) 


This  expression  is  plotted  in  Fig.  3  for  different  pressures  (different 

values  of  n  ).  The  increase  and  eventual  decrease  with  increasing  H., 

8 

pressure  is  clearly  seen  in  the  figure. 

For  a  fixed  x  in  Fig.  3,  the  precursor  electron  number  density  as  c. 
function  of  partial  pressure  can  be  obtained.  This  is  plotted  in  Fig.  4 
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for  x  =  50  cm.  and  x  =  cm.  These  curves  show  the  "peak  behavior"  and  the 
partial  pressures  at  which  the  peaks  occur  agree  with  the  values  in  Fig.  2, 
which  were  obtained  from  a  more  general  analysis. 

In  Fig.  5,  a  comparison  is  made  between  the  precursor  profile 
obtained  from  this  modified  theory  (at  a  partial  pressure  of  50|iHg.  ) 
and  the  original  Wetzel  theory^  (assuming  a  one-step  photoionization  of 
H^).  The  electron  number  densities  obtained  from  the  original  Wetzel 
theory  are  several  orders  of  magnitude  higher  and  the  fall-off  with  distance 
ahead  of  the  shock  much  less  rapid  than  obtained  from  the  modified  theory. 
The  experiments  n  Ref.  18  agree  quite  well  with  the  modified  theory  but 
not  at  all  with  the  original  Wetzel  theory. 

If  In  n(xj  is  obtained  from  Er.  39.  and  if  this  expression  is  then 
differentiated,  the  following  equation  for  the  slope  of  the  logarithm  of  the 
precursor  profile  results: 


in  n(x) 


x 


1  + 


I  +  V 


16.  3 


[  1  +  e  '  t|  Ei(  -  r|  *]  +-|^3 


(40) 


Fig.  6  is  a  plot  of  this  expression  as  a  function  of  the  H,  partial  pressure 

C* 

for  a  fixed  distance,  x  =  100  cm,  ahead  of  the  shock.  This  curve  will 
be  compared  with  experiment  in  Ref.  18. 
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SECTION  V 


CONCLUDING  REMARKS 

Until  recently,  very  little  experimental  data  on  electron  precu~sors 
v  18  available.  The  recent  shock  tube  experiments  in  Argon^’  *3  provided 
data  on  which  to  test  the  several  theories  which  have  been  introduced  over  the 
past  few  years.  For  one  thing,  the  experiments  of  Lcderman  and  Wilson 
strongly  suggest  that  "diffusion  theories",  at  least  for  the  their  experiments; 
do  not  apply.  They  showed,  experimentally,  that  radiation  from  the  hot 
gas  behind  the  shock  front  was  responsible  for  the  observed  ionization  ahead 
of  the  shock.  They  also  found  an  electron  density  distribution  ahead  of  the 
shock  which  did  not  change  with  time  in  a  shock  fixed  coordinate  system. 
(Holmes*3  did  not  fmd  this  "steady  state"  situation  in  his  experiments.) 

These  experiments  strongly  recommend  a  photoionisation  theory  like 
the  one  proposed  by  Wetzel3.  Wetzel' s  theory  postulates  a  plane 
source,  where  in  most  applications  (including  the  present  shock  tube  application) 
a  point  or  spherical  radiator  is  more  realistic.  Moreover,  when  the  precursor 
number  densities  due  to  photoionization  of  the  Argon  driven  gas  is  calculated 
by  Wetzel' s  theory,  the  levels  are  found  to  be  several  orders  of  magnitude 
below  the  experimental  values. 

With  these  limitations  of  the  Wetzel  theory  in  mind,  the  present 
modifications  of  his  theory  were  rrade: 

1)  A  point  or  spherical  radiator,  located  at  the  shock  front  and 
moving  with  the  shcck  velocity,  was  assumed. 

2)  The  ionization  ahead  of  th^  shock  was  viewed  not  as  a  one- step 
photoionization  of  Argon  but  as  a  one-step  photoionization  of  impurities 
present  in  the  Argon  driven  gas. 

3)  Since  the  radiation,  for  frequencies  below  the  f-'rst  ionization  edge 


22 


of  Argon,  w.as  not  necessarily  expected  *o  behave  like  black  boor,  radiation, 
this  question  was  investigated.  It  was  found  that,  because  of  the  impurities 
in  the  driven  gas,  a  black  body  radiator  is  a  reasonable  one  to  assume. 

As  an  application  of  the  theory,  precursor  profiles  were  calculated 

for  a  one -step  photoionization  of  an  H.,  impurity  (the  partial  pressure  of 

was  varied  from  2  to  75  uHg.  ),  The  driven  gas  was  assumed  to  be 

Argon  at  a  pressure  of  2  mm.  Hg.  ;  and  the  Mach  number  was  assumed  to 

be  13.  2,  giving  a  temperature  of  approximately  11,  0G0°K  behind  the  shock. 

A  comparison  is  made  with  the  original  Wetzel  theory,  ap  Alied  to  a  one- 

step  photoionization  of  H^.  The  electi  on  number  densities  obtained  from 

the  Wetzel  theory  are  several  orders  of  magnitude  higher  and  the  fall-off 

with  distance  ahead  of  the  shock  much  less  rapid  than  obtained  from  the 

modified  theory  (Fig.  5).  Also,  according  to  the  Wetzel  theory,  the 

precursor  electron  number  density  would  decrease  with  an  increase  in 

partial  pressure,  whereas  the  modified  theory  predicts  an  increase 

and  eventual  decrease  with  increasing  partial  pressure.  Experiments 

have  been  performed  where  minute  quantities  cf  were  added  to  the  Argon 

18 

driven  gas  and  changes  in  the  precursor  profile  were  recorded 

In  all  respects,  the  experiments  agree  quite  well  with  the  modified  theory 

but  not  at  all  with  the  original  Wetzel  theoiy. 
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